UNCLASSIFIED 

m  410213 

DEFENSE  DOCUMENTATION  CENTER 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 

CAMERON  STATION.  ALEXANDRIA.  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  vlth  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other- 
vise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  Invention  that  may  in  any  way  be  related 
thereto. 


*,  ■  fr  a  z 


MEMORANDUM  REPORT  NO.  14<i 
APRIL  1963 


ON  THE  QUASI-LINEAR  SUBSTITUTION  METHOD  FOR 
MISSILE  MOTION  CAUSED  BY 
STRONGLY  NONLINEAR  STATIC  MOMENT 

Charles  H.  Murphy 


RDT  &  E  Rraiacl  No.  1M010501A005 

I  BALLISTIC  RESEARCH  LABORATORIES 

ABERDEEN  PROVING  GROUND,  MARYLAND 


ASTIA  AVAILABILITY  NOTICE 

Qualified  requestors  nay  obtain  copies  of  this  report  from  ASTIA. 


The  findings  In  this  report  are  not  to  be  construed 
as  an  official  Department  of  the  Any  position. 


BALLISTIC  RESEARCH  LABORATORIES 


MEMORANDUM  REPORT  NO.  1466 


APRIL  1963 


ON  THE  QUASI- LINEAR  SUBSTITUTION  METHOD 
FOR  MISSILE  MOTION  CAUSED  BY 
STRONGLY  NONLINEAR  STATIC  MOMENT 


Charles  H.  Murphy 


Exterior  Ballistics  Laboratory 


RDT  &  E  Project  No.  1M010501A005 


ABERDEEN  PROVING  GROUND, 


MARYLAND 


BALLISTIC  RESEARCH  LABORATORIES 


MEMORANDUM  REPORT  NO.  1466 


CHMurphy/ Jdk 

Aberdeen  Proving  Ground,  Md. 
April  1963 


ON  THE  QUASI- LINEAR  SUBSTITUTION  METHOD 
FOR  MISSILE  MOTION  CAUSED  BY 
STRONGLY  NONLINEAR  STATIC  MOMENT 

ABSTRACT 

An  improved  quasi-linear  substitution  method  is  developed  to  treat  properly 
the  influence  of  a  cubic  static  moment  on  the  modal  damping  of  a  missile  acted 
on  by  quite  general  nonlinear  damping  and  Magnus  moments.  The  predictions  of 
this  method  are  compared  for  various  special  cases  with  those  of  the  more 
accurate  but  much  more  complicated  perturbation  method.  The  new  quasi-linear 
theory  predicts  boundary  curves  for  planar  motion,  almost  circular  motion  and 
almost  planar  motion  which  are  quite  close  to  those  of  the  perturbation  theory. 
An  original  result  of  the  theory  is  that  all  planar  singular  points  for  a  non¬ 
spinning  missile  whose  moment  coefficients  are  only  functions  of  the  total  angle 
of  attack  are  nodes.  That  is,  almost  planar  motion  with  amplitude  close  to  that 
of  a  stable  planar  limit  motion  will  tend  to  that  motion. 
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Subscripts 

c  quantities  evaluated  for  circular  motion 

p  quantities  evaluated  for  planar  motion 
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1.  INTRODUCTION 


In  Reference  1,  three  different  quasi-linear  methods  were  described  and 

their  different  predictions  of  the  nutational  frequency  for  a  missile  with  a 

cubic  static  moment  were  compared  with  the  exact  result  obtained  by  the  use 

of  elliptic  integrals.  The  best  of  the  three  was  called  the  substitution 

method  and  was  employed  to  obtain  the  combined  effect  of  a  cubic  static  moment , 

varying  air  density  and  both  linear  and  nonlinear  damping  moments.  Although 

2  3 

this  approach  was  not  as  accurate  as  the  perturbation  method  which  uses  the 
exact  elliptic  function  solution  for  no  damping  as  the  initial  approximation, 
it  did  give  trends  with  a  significant  reduction  in  the  necessary  algebraic 
work. 


One  difficulty  with  the  substitution  method  was  use  of  a  rather  strange 
condition  on  the  damping  of  the  modal  amplitudes: 


,  .  i02 

K  e  +  I^e  *  0 


(1) 


A  re-examination  of  this  question  has  revealed  a  different  substitution,  method 
which  yields  the  same  expression  for  the  frequencies  but  does  not  make  use  of 
Equation  (l).  This  improved  substitution  method  does  lead  to  different  ex¬ 
pressions  for  modal  damping  which  are  the  same  as  those  for  the  earlier  substi¬ 
tution  method  for  planar  motion  but  are  closer  to  those  obtained  from  the 
perturbation  method  for  almost  circular  motion.  It  is  the  purpose  of  this 
report  to  describe  this  new  quasi-linear  method. 

2.  THE  IMPROVED  SUBSTITUTION  QUASI-LINEAR  SOLUTION 

The  equation  for  the  pitching  and  yawing  motion  of  symmetric  missiles  can 
2 

be  written  in  the  form 

1"  +  (H  -  h  -  iP)T'  -  (M  +  iPT)I  =  0  (2) 


where  the  various  symbols  are  defined  in  the  Table  of  Symbols . 
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(?) 


f  +  (Ho  -  iP)T’  -  (Mo  +  iPTo)I 

*  -  [h  -  H0  -  {'  +  £(M  -  Ho)  +  iP(T  -  Tq)J  T 

where  the  zero  subscript  denotes  the  value  of  the  aerodynamic  quantities  for 
zero  amplitude  motion. 


The  left  side  of  Equation  (3)  is  the  linearized  version  of  Equation  (2)  and 
has  a  solution  in  the  form 


i02 

|  =  Kxe  A  +  Kge 


(k) 


The  primary  interest  of  a  missile  designer  is  usually  the  behavior  of  the 

modal  amplitudes  K..  This  can  conveniently  be  described  by  their  logarithmic 
J 

derivatives 


t 


(5) 


For  linear  moments,  the  X*s  are  constants  but  when  the  moments  are  nonlinear 

u 

the  various  quasi-linear  analyses  obtain  which  are  functions  of  and 

K„. 


We  now  differentiate  Equations  (4)  twice,  substitute  in  Equation  (3)  and 
solve  for  the  frequency  and  damping  of  the  first  mode. 


<^>2  +  *!<*!  +  V  -  h 

-  i  [(20'  -  P)X1  +  -  PT0  *  0"] 

■[H'E0  -f]  [H  *  4>  ♦  (»a  ♦  i^)  ^  «-1?] 

-  f(M  -  M0)  +  iP(T  -  ijj  [l  +  ^  e-^l  (6) 
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(6) 


-{[■ 


(0o)2  -  P0O  +  M 


-  i 


”(202 


P)\c 


where 


i  ■  0,  -  A 


For  linearized  motion,  all  the  terms  on  the  right  of  Equation  (6)  vanish 
except  for  the  term  in  braces.  Since  the  left  side  of  the  linearized  Eq.  (6) 
is  constant  or  slowly  varying  due  to  density  or  Mach  number  induced  variations 


of  M  , 
o' 


Hq,  and  Tq  and  spin  or  velocity  variations  in  P,  it  can  only  equal  the 


1? 


periodic  term  in  e  r  on  the  right  side  if  both  sides  are  zero.  This  condition 
yields  the  usual  linear  relations.  The  nonlinear  terms  on  the  right,  however, 
will  in  general  be  periodic  but  their  average  or  d.  c.  components  will  not 
necessarily  be  zero.  This  average  will  directly  affect  the  terms  on  the  left 
side  and  it  is  the  assumption  of  the  quasi-linear  method  that  the  average  is 
the  only  influence  of  the  nonlinearities  on  the  frequencies  and  damping 
exponents.  We,  therefore,  average  Equation  (6)  over  a  period  of  nutation*  and 
neglect  the  small  damping  term  in  comparison  with  Mq  in  the  real  part  of  this 
equation, 

(0^)2  -  P0[  +  Mq  -  i  -  P)X1  +  -  PTo  +  0"J 

-  k  f*  <H  -  %  -  f >  [<  +  4  (jj)  *'1?]  * 

■  h  £  [<M  -  M»>  ♦  iP<T  -  v]  [x  *  \  «‘1?]  * 


7  which  is  the  cosine  of  the  total  angle  of  attack  can  be  related  to 
|  |  j  ,  the  magnitude  of  the  sine  of  the  total  angle  of  attack  by  the  relation 


Ballisticians  frequently  use  the  terms  nutation  and  precession  to  distinguish 
two  modal  oscillations.  Nutation  in  this  report  has  the  classical  meaning 
assigned  by  top  theory,  i.e.,  the  variation  of  the  amplitude  of  the  total 
angle,  |  5 |  • 
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(8) 


2  c2 
y  *  1  -  8 


where  62  **  |  i  j  2  ■  5  T 


6  can,  then,  be  computed  from  Equation  (k) 
52  *  +  K^Ce1^  +  e'if^) 

=  +  ^Kg  cos  0 


Since 


21 

7 


.  ill  m 

1>7  2(1  -  o  ) 


(9) 


(10) 


we  see  that  y  fy  is  an  odd  function  of  $  and,  therefore,  only  affects  the  real 
part  of  Equation  (7).  If  we  make  the  assumption  that  H  and  T  are  functions  of 
6  2  and,  thereby,  are  even  functions  of  %,  the  following  equations  may  be  obtained 
from  Equation  (7). 


Similar  expressions  apply  for  the  other  mode. 


Ik 


This  result  differs  from  that  of  the  substitution  method  of  Reference  1 

j 1  * 

in  the  presence  of  20  -  P  in  the  denominators  of  X,  and  X,  .  In  that  report, 
ii  X  i  x  x 

01  -  0g  appeared.  This  is  the  same  as  20^  -  P  for  planar  motion  but  differs 
for  other  motions. 

i 

For  small  geometrical  angles  (7  *  0)  and  a  cubic  static  moment 

2 

(M  =  Mq  +  Mg6  ),  Equation  (ll)  and  the  corresponding  equation  for  the  other 
mode  yield  frequency  relations  which  are  identical  with  those  of  Reference  1. 


Equation  (13)  can  he  differentiated  and  substituted  in  Equation  (12) 
with  the  result: 


(15) 
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A  similar  equation  can  be  derived  for  the  other  mode. 


(16) 

Equations  (15-16)  may  now  be  solved  simultaneously  for  \ ^  and  \r>. 


(17) 

where  the  a.  ,'s  are  defined  in  the  Table. 
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TABLE  (Cont.) 


/ 
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5.  PLANAR  MOTION 


In  the  Table,  values  of  the  a.,  's  are  computed  for  quasi-planar  motion, 

i.e.,  nutation  between  zero  and  b  *  For  this  motion,  K.  =  =  K,  & 

max  9  1  z  7  max 

»  +  Kg  a  2K  and  the  ratio  of  the  maximum  value  of  the  cubic  term  to  the 

linear  term  is 


m  = 


yL 

Vmax 


(18) 


There  are  three  types  of  cubic  static  moments  which  can  cause  a  periodic 

2 

motion*  These  are  shown  in  Figure  1  and  may  be  described  as  follows  : 

(a)  Stable  at  small  angles;  more  stable  at  larger  angles  <  0,  Mg  <  0) 
quasi-planar  motion  of  all  amplitudes  possible  (m  >  0) . 

(b)  Stable  at  small  angles;  less  stable  at  large  angles  (Mq  <  0,  Mg  >  0) 
amplitude  of  quasi-planar  motion  must  be  small  enough  to  assure  positive 

M(  -  1  <  m  <  0)  • 

(c)  Unstable  at  small  angles;  stable  at  large  angles  (ft  >  0,  Mg  <0)J 
amplitude  of  quasi-planar  motion  must  be  large  enough  to  ensure  positive 
average  M  (m  <  -  2). 

Pure  planar  motion  occurs  when  spin  is  zero.  With  zero  spin,  a  number 
of  simplifications  are  possible. 
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62  »  2K2(l  +  cos  % 

(25) 

(82)'  =  -  ^K2  sin  % 

(214) 

h  a  X2 

2(14  +  3m)X*  -  8  ^  )  -  6m  ^  ) 

8  +  9m 

(25) 

If  the  cause  of  the  variation  in  the  coefficients  is  changing  air  density  due 
to  entering  or  leaving  an  exponential  atmosphere^, 


i  %  £ 

M0  "  ^  ”  P 


a 


•  •  —  Xg  — 


*(*♦*)  fx*-f) 


8  +  901 


(26) 


(27) 


/  X  k 

Equation  (27)  was  essentially  derived  by  Coakley,  Laitone  and  Mass  and 

predicts  that  the  amplitude  of  planar  motion  for  a  type  (b)  moment  has  an 

upper  bound  imposed  by  the  requirement  -  8/9  <  m.  For  linear  damping,  however, 

the  more  exact  perturbation  method  of  Reference  2  yields  the  relation:1'^ 


-  bQ 

X1  =  X2  =  2(l  +  m) 


(28) 


where 

bQ  =  4(1  +  m)a2  -  2ma^ 


types  (a)  and  (c) 


■  2(2  +  m  -  2a2  -  ma^) 


type  (b) 
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S2  *  ‘p2 


E1  -  vkp] 


“j,  ■  <V5>  *'2  [20  +  kp)»2  -  El 


K  =  K(k  )  complete  elliptic  intej^al  qf  f^pst  kind 


E  =  E(k  )  complete  elliptic  integral  of  the  second  kind 
P  P 


.2 _ m 

kp  "  2(1  +  m) 


-m 

~  2  +  m 


types  (a)  and  (c) 
type  (b) 


[M] 


The  coefficients  of  l^  +  il  in  Equations  (27)  and  (28)  are  compared 
for  the  three  types  of  moments  in  Figures  (2-4).  With  the  exception  of  the 
vicinity  of  m  ■  -  8/9,  the  quasi-linear  substitution  value  is  a  reasonable 
approximation  of  the  more  exact  but  quite  complicated  perturbation  result. 

4.  ALMOST  CIRCULAR  MOTION 

For  almost  circular  motion,  one  modal  amplitude  is  much  larger  than  the 
other  and  any  static  moment  can  be  approximated  by  a  cubic  in  the  vicinity  of 
the  amplitude  of  the  circular  motion.  In  the  Table,  the  coefficients,  a.,  , 
are  computed  for  Kg  <  <  K^.  (The  coefficients  for  <  <  Kg  may  be  obtained 
by  interchanging  and  Kg.)  As  in  the  case  of  planar  motion,  a  number  of 
simple  relations  can  be  written  for  Kg  <  <  and  constant  spin.*  (P  »  0) 


-  |  +  V  -  ftQ(l  +  m) 


(29) 


1  _  JZ 

2  V 


(1  +  2m) 


(50) 


If  KL  <  <  Kg 


,  the  frequency  equations  are 

*  |  +  J  -  Mo(l  +  2m);  $2~\~ 


/ 


7T. - - 

-  Mo(l  +  m) 
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2  it 


Note  that  Equation  (32)  contains  the  very  large  quantity  K^/K^.  For  most 

nonlinearities ,  the  averaging  process  of  the  integral  formally  yields  a 
*  -2 

magnitude  for  of  the  order  of  10  and  certainly  much  less  than  one.  This 
need  for  the  average  of  a  large  periodic  term  to  vanish  would  lead  us  to  .expect 
Equation  (32)  to  be  less  accurate  than  Equation  (31) •  As  we  shall  see,  this 
is  a  correct  conjecture. 
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on  m: 


The  exact 


elliptic  integral  solution 


2 


places  the  following  limitations 


type  (a)  no  limitation  (m  >  0). 

type  (b)  only  possible  circular  motions  are  those  for  which 
-  2/3  <  m  <  0  . 

type  (c)  only  possible  circular  motions  are  those  for  which 
M  is  negative  (m  <  -  l). 

The  presence  of  2  +  3m  in  the  denominator  of  Equation  (33)  is  the  first  time 
a  quasi-linear  approach  has  been  able  to  indicate  the  completely  unexpected 
limitation  on  circular  motions  for  a  type  (b)  moment  which  was  previously 
obtained  through  the  use  of  elliptic  integrals.  This  result  is  the  first 
evidence  of  the  value  of  the  improved  substitution  method  in  comparison  with 
that  of  Reference  1.  The  fallacious  indication  of  trouble  for  m  =  -  l/2  in 
Equation  (3*0  reinforces  our  concern  *or  the  value  of  the  expressions  for 
damping  of  the  small  modal  amplitude. 

In  order  to  derive  an  estimate  for  the  accuracy  of  Equations  (33-3^) >  we 
will  now  consider  two  special  cases  for  which  the  coefficients  are  constants 
(m’  =  Mg  =  0). 


X 


1 


2  +  2m  ,* 

5  +  3m  X1 


(35) 


*  _  -  2m(l  +  m) 

X2  *  (2'+"3m)  (1  +  2m) 

These  cases  were  treated  by  the  perturbation  method  in  Reference  2  and  the 
validity  of  our  results  will  be  determined  by  comparison  with  results  of  that 
method « 


^1  +  Xg 


(36) 
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4.1  CUBIC  MAGNUS  MOMEWT 


In  the  first  example,  we  will  consider  a  spinning  missile  with  a  cubic 

2 

Magnus  moment  (T  *  Tq  +  TgB  ) 


*  -  H 

*  o 


X1  = 


[j  +  7-  Mq(1  +  m)J  +  P  [to  +  TgK^I 
2  *f-  A  (l  +  m) 


■  -[l  '2Ho  -  |1+"> 


-V2  A 

P(2T  -  H  +  2T, 
x  o  o  t 


*1 


(37) 


H.  I 


X2  ~ 


y-6(i  +  a,)3- 


2  */-  M  (1  +  2m) 


2H  +  1  +  2m 

o 


-1/2 


P(2T  -  H  +  4T0, 
oo  2 


A  p 

where  P  » 


N 


175 


(38) 


When  the  damping  coefficient  of  the  large  mode  as  given  by  Equations  (35) 
and  (37)  is  compared  with  that  obtained  from  the  perturbation  method  ,  we  find 
them  to  be  identical l  A  comparison  of  the  damping  coefficients  for  the  small 
mode  reveals  that  they  do  differ.  A  measure  of  the  magnitude  of  this  difference 
can  be  obtained  by  considering  the  conditions  for  a  circular  motion  singularity 
which  is  a  stable  node.  The  location  of  the  singularity  is  given  by 

4  ■  -(!)(Rte)  [2Ho  -  I1  +  "  |  1/2  H  -  Ho  +  ■ 0 

(39) 


24 


Since  small  circular  motion  must  grow  and  large  circular  motion  must  decay, 


2H  -  1  +  m 
o 


pt2<0 


-1/2  A 

P(2Tq  -  Hq)  <  0 


(40) 

(41) 


These  inequalities  naturally  are  equivalent  with  those  of  Reference  2,  The 
final  requirement  is  that  in  the  vicinity  of  the  singular  point,  almost 
circular  motion  will  approach  circular  motion. 

\2  =  \*2  <  0  (42) 
By  the  use  of  Equations  ( 38-39) >  Tg  may  be  eliminated  from  Inequality  (42). 

2  [lt2  rfal  l/2J  Ho-|i  +  a»|'1/2$<2To-Ho)>0 


(43) 


But  Inequalities  (40)  and  (43)  for  m  outside  the  interval  (-1,  -l/2)  may  be 
combined  to  yield 

1/2  A  f.  .1/2  ,  ,1/2~1 


2  1  +  m  H  <  P(2T  -  H  )  <  2 

o  o  o 


|l+2m  +  2  1  +  m  J J  Hq 


(44) 

Outside  the  forbidden  interval  of  -  1  <  m  <  -  l/2  Inequality  (44)  requires  that 

1/2 


Ho  be  positive  and 


,  l/2  P(2T  -  H  ) 

<— i-2 


-L/* 

1  +  2m  +2 

1  +  m 

(45) 
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The  corresponding  inequality  from  Reference  2  is* 


1  +  m 


1/2 


s, 

P(2T 


o 

SF 


Ho> 


16 


7m  | 


1  +  m 


-1/2 


•  o 

These  bounds  sure  compared  in  Figure  5*  As  can  be  seen  from  the  figure,  the 
upper  bound  is  a  good  approximation. 


(46) 


4.2  CUBIC  DAMPING  MOMENTS  AND  ZERO  SPIN 

5 

For  a  nonspinning  missile,  it  has  been  shown  that  there  are  two  cubic 
damping  moments  which  can  affect  the  modal  amplitudes:** 


H  »  Hq  +  (Hg  +  MU)6C 


M  =  Mq  +  M 2&2  +  Mn(62)' 


(*7) 

(U8) 


Since  we  are  concerned  with  almost  circular  motion  (K^  >  >  Kg)>  a  much  more 
general  moment  can  be  considered.  This  moment  will  be  approximated  by 
Equations  (47-48)  for  almost  circular  motion.  This  moment  can  be  written  in 
terms  of  the  aerodynamic  moment  coefficients  as*** 


C 

m 


+  iC 

n 


i 


+  c2  6 


+  c*)  i  +  d 


(49) 
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1/2 


The^  of  that  report  is  |  1  +  m  |  P  for  m  >  -  2/3  and  -  P  1 1  +  m  | 
for  m  <  -  1. 

ii  || 

The  presence  of  M.  in  (1*7)  is  due  to  its  definition  as  the  coefficient  of 
2  ,  -L-l 

|  I  in  Reference  5* 

The  good  approximation  i  =  i  (q  +  ir)  iV"  has  been  made  in  Equation  (49) 

so  that  Cw  and  C.,  appear  as  a  sum. 

M  M. 

q  a 
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*  */  2  2  1 
where  c  =  c  f(fi  )  J  is  a  function  of  (6  ) 

d  *  d  (6^)  is  a  function  of  6^ 

*.  . 

c  (0)  =  0  and 


a  (o)  -  cM  + 

10 


The  coefficients  of  the  differential  equation  of  pitching  and  yawing 
motion  assume  the  form 


+  C2® 


*(  <^>) 


(50) 

(51) 


If  H  and  M  are  differentiable  functions,  they  can  be  expanded  about  the 

2  1 

circular  motion  with  amplitude  &c  and  amplitude  derivative  (6  )  «  0. 

(52) 

(53) 

where  H  -  H(B^) 
c  c 
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62  -  +  2KX Kg  cos  $ 

(B2)'  -  -  -  0g)  8i“^ 


Equations  (52-53)  are  essentially  of  the  same  form  as  Equations  (1*7-48)  but 
allow  us  to  consider  much  more  complicated  moments. 


For  these  moments  and  almost  circular  motion, 


(55) 


The  actual  damping  exponent  for  the  larger  mode  can  be  obtained  from  Equations 
(35)  and  (54)  and  is  exactly  that  given  by  the  perturbation  method. 


(56) 


The  conditions  on  for  a  stable  node  at  K.  ®  6  are 
l  1  c 

Hc  -  O  (57) 

Hc>°  <*> 


According  to  Equation  (36)  the  actual  damping  of  the  smaller  mode  near  the 
*  ,  * 

singularity  (x.  ■  0)  Is  X2*  For  a  stable  node  this  must  be  negative. 
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The  corresponding  inequality  derived  from  the  perturbation  method  is 


These  upper  bounds  are  compared  in  Figure  6.  It  is  interesting  to  note  that 
#  # 
the  M  function  is  necessary  for  a  circular  limit  cycle,  i.e.  M  *  0  does 

not  satisfy  Inequality  (6l). 


5.  ALMOST  PLANAR  MOTION 

Another  important  special  type  of  motion  is  almost  planar  motion  (K^  •  Kg)* 
For  this  case  somewhat  more  lengthy  algebra  is  required  but  fairly  general 
results  are  attainable.  For  planar  motion  and 


62  =  21^(1  +  cos  %) 

P  P 

(62) 

sl“  $ 

(63) 

0  •  J-  M  (4  +  3m  ) 

o  p 

(64) 

„ . 

where  m  =  — zt-*- 

* 

If  H  and  M  are  differentiable  functions,  they  can  be  expanded  about  planar 
motion  with  amplitude  2K^. 
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(65) 


(66) 


where  the  p  subscript  denotes  quantities  evaluated  for  planar  motion  with 

O  P  ,  P  I  PI 

6  =  8p  and  (8  )  «  (5  )p  . 


For  a  nonspinning  missile, 


(67) 


In  order  to  determine  the  character  of  a  planar  singularity  amplitude  plane, 
differential  equations  should  he  derived  for  the  neighborhood  of  the  singularity. 
The  variables  and  are  introduced  for  this  purpose  and  squares  of  these 
variables  will  be  omitted  in  comparison  with  the  variables  themselves. 


K>  =  Kp(l  +  6j) 


82  -  6p  =  21^(1  +  cos  ?)  (Cl  +  e2) 


(68) 

(69) 


%'  *  -  02 


.(1/2)  [■/-  Mo  [4  +  mp(5  +  2^  +  1|€2)  J  +J-  M0  [j*  +  mp(5  +  4^  +  2e2)]  I 


3'  fi* 

p  L 


(62)’  -  (62)’ 


3m  “I 

«  +  6m  ^€1  +  €2M 

P  —* 

sin  0  J  -  jj  2K^p  sin 


(70) 


K2^’  sin  % 
P  P  ^ 


f8t9°A  1 

{rrjrJ  <ei 


+  e2) 


(71) 
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result  reduced  by  routine  algebra  so  that  Equation  (75)  assumes  the  much 
simpler  form 


31 


* f  [*], 


cos 


J  -  0,2 


2n 

Mi  -I  f 

J  O 


H  cos 
P 


Similarly 


*  * 

Xrt  *  X_  m 


2  ■  *P  "  Ml.  'rjry  {[Vi  *  <8  +  +  [(l6  +  “V*i  *  r< 

-  [TV1  +  (l6  +  17"p)s2]  r2  |  (75) 


The  numerical  subscript  on  the  outside  of  the  bracketed  expressions  in  the 

definition  of  X.  identifies  that  expression  as  a  particular  Fourier  cosine 
J 

coefficient.  It  is  quite  surprising  that  the  influence  of  H  on  X  is  completely 

1  * 

determined  by  its  first  order  Fourier  cosine  coefficient.  The  influence  of  M  , 

however,  is  specified  by  the  zeroth  and  second  order  Fourier  cosine  coefficients 

of  its  first  derivative.  These  coefficients  are  computed  for  fixed  modal 

amplitudes,  K.,  and  therefore,  are  functions  of  these  amplitudes. 

J 

* 

For  a  planar  singularity  Xp  ®  0  and  Xj  can  be  computed  from  the  following 
special  form  of  Equation  (17) 


ajlXl  +  SJ2X2 


(76) 


where  a, 


11 


*22 


2(4  +  3m  )(8  +  7m) 


al2  •  a21 


-  4m  (4  +  3m ) 

W+frJVS  +  5mJ 
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(77) 


DQ 


ll  ■  ■  5flTP  J^r«r+Ty  [bsi  *  “2] 


M 


l2  *  •  4(8  ;'9^kb  v'sji  0*1 +  ”'2] 


(78) 


2 

where  a=-m(8+3m)+  (128  +  200m  +  75m  )  r 

P  P  P  p  0 

-  (24mp  +  15mp)  r2 

b  =  64  +  96mp  +  55mp  -  3mp(8  +  5mp)  rQ 

-  (128  +  248mp  +  105znp)  r2 

The  differential  equation  for  solution  curves  in  the  vicinity  of  a  planar 
singularity  in  the  amplitude  plane  is 


de2  ^2  ael  +  l362 

d£i  X1  ”  be-^  +  ac2 


(79) 


According  to  Reference  6,  the  singularity  must  be  either  a  saddle  point  or  a 

2  2  2  2 
node.  It  is  a  node  if  a  -  b  is  negative  and  a  saddle  if  a  -  b  is  positive. 

a^-b^=  (a  +  b)(a-b) 

=  -  4(8  +  5mp)(8  +  9mp)  £2(2  +  mp)  -  (rQ  +  r2)(8  +  5mp)j 
x  [(4  +  3mp)(l  +  2rQ)  -  4(2  +  3mp)  r2Q  (80) 


.Note  that  if  r  and  r0  vanish  and  m  is  outside  the  interval  (-2,  -  8/9), 

2  2  o  d  p 

a  -  b  is  negative.  Thus,  if  the  aerodynamic  moment  coefficients  are 
2  2  1 

functions  of  6  alone  and  not  functions  of  (6  )  ,  all  planar  singularities 
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are  nodes  and  almost  planar  motions  near  a  planar  singular  motion  will  tend 
to  the  planar  singular  motion  If  neighboring  planar  motions  tend  to  the  planar 
singular  motion.* 

Another  interesting  special  case  is  that  when  the  moment  coefficients 
are  precisely  those  defined  by  Equations  (1*7-48) 


H1  -2(H8*«u)^ 


(81) 

(82) 


(83) 

(84) 

(85) 


Therefore,  a  planar  singularity  is  a  node  if 


(8  +  9mp)(8  +  5mp)(4  +  3mp)(l  +  2rQ)  £l*  +  2mp  -  rQ(8  +  5mp)J  >  0 


(86) 


When  m 


to 


is  outside  the  interval  (  -  2,  -  8/9)  this  Inequality  is  equivalent 


-  1  < 


-  M, ,  1*  +  2m 

<  ,  P 

Hg  *  15  +  7mp 


(87) 


For  circular  singularities,  Equation  (60)  shows  the  circular  singularities 

2 

are  always  saddles  if  the  moment  coefficients  are  functions  of  6  alone, 


i  •  e  • , 


“u 


dM 

di?j 


* 


c 


3k 


In  Reference  3,  the  perturbation  method  was  applied  to  almost  planar  motion 
and  after  considerable  tedious  algebra,  inequalities  like  (87)  were  obtained. 
The  lower  bounds  are  identical  but  the  perturbation  method* s  upper  bound  is 
expressed  in  terms  of  complete  elliptic  integrals  and  differs  from  that  of 
Inequality  (87)  when  m^  ^  0.  These  two  boundary  curves  are  compared  in 
Figure  7.  The  much  more  easily  derived  bound  of  Inequality  (87)  is  surpris¬ 
ing!  y  good. 


SUMMARY 


1.  A  quasi-linear  substitution  method  has  been  derived  on  plausible 
assumptions  and  compared  with  the  more  exact  results  of  the  more  complicated 
perturbation  method. 

2.  The  planar  motion  predictions  are  good  when  m  is  not  near  -  8/9, 

3.  The  damping  of  the  dominant  mode  of  almost  circular  motion  is  exactly 
predicted  while  the  damping  of  the  other  mode  yields  approximately  correct 
stability  boundaries. 

4.  The  character  of  planar  singularities  can  be  reasonably  well  deter¬ 
mined  by  this  method. 

5.  In  view  of  the  above,  the  algebraically  much  simpler  quasi-linear 
method  can  be  used  to  obtain  approximate  stability  boundaries  in  the  presence 
of  a  strongly  nonlinear  static  moment. 
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FIG.  2 


COMPARISON  OF  PLANAR  DAMPING 
MOMENT  COEFFICIENTS 
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